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Abstract. We identify a duality transformation in one-dimensional hopping 
models that relates propagators in general disordered potentials linked by an up- 
down inversion of the energy landscape. This significantly generalises previous 
results for a duality between trap and barrier models. We use the resulting insights 
into the symmetries of these models to develop a real-space renormalisation 
scheme that can be implemented computationally and allows rather accurate 
prediction of propagation in these models. We also discuss the relation of this 
renormalisation scheme to earlier analytical treatments. 
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1. Introduction 

The motion of particles in disordered environments is important in many contexts, 
from glass- forming liquids and colloids [H j to biomolecules moving in the crowded 
environment of the cell [3], to electrical properties of disordered materials 4J. In this 
article, we discuss subdiffusive propagation in simple one-dimensional models. While 
the case of one-dimensional motion may seem simplistic, it is relevant for a variety of 
model systems: from early studies of electrical transport |4J to recently-defined models 
of glassy behaviour [S], and also to the motion of defects in disordered magnets, to 
disordered elastic chains and to networks of resistors and capacitors (see [51 [71 13 IS] 
for reviews). 

The results that we will present are based around a duality symmetry which 
relates pairs of discrete one-dimensional models. In a previous study [10] , we showed 
that motion in apparently disparate models can be related exactly, at fixed disorder. 
Here we generalise those results to a much wider range of hopping models, by treating 
their master equations in a simple operator formalism. We also discuss how these 
results are related to earlier studies of disordered reaction-diffusion systems by Schiitz 
and Mussawisade [TT] . Based on the symmetries of the problem, we then introduce 
a real-space renormalisation scheme in the spirit of that of Ic Doussal, Monthus and 
Fisher T2]: our scheme is implemented computationally and allows rapid prediction 
of the propagation in these energy landscapes for fixed disorder, at low computational 
cost. 

The central feature of the renormalisation scheme is that on a given time scale, we 
have a procedure for decomposing the system into effective trap and barrier regions. 
Particles within trap regions equilibrate there, while those in barrier regions decay 
into the effective traps. The duality symmetry relates trap and barrier regions of 
pairs of models, and demands that they be treated on an equal footing within the 
renormalisation scheme. In some sense, the renormalisation scheme is connected with 
the ideas of an energy landscape in these disordered systems [I^ , but we note that all 
properties of the energy landscape are here derived directly from the master operator 
of the stochastic dynamics. In this sense, our scheme allows the energy landscape 
to be derived from the dynamical rules of the system: this is the opposite of the 
usual situation in which thermodynamic properties are used to infer the routes by 
which dynamical processes take place. Thus, while our results are clearly restricted 
to a very simple class of models, it is natural to ask if they might be generalised to 
higher-dimensional energy landscapes. 

The form of the paper is as follows: in Sec. [21 we define our models and give 
the duality relation between their master operators. The consequences of the duality 
relation for propagation in these models are discussed in Sec. [S] In Sec. [H we explain 
our effective dynamics scheme; Sec. [5] contains numerical results for specific ensembles 
of disordered models; and Sec.[6lcloses with a brief summary and some open questions. 

2. Models and duality relations 

We define a disordered one-dimensional hopping model in terms of rates for hops from 
site i to sites i — \ and i -I- 1, which we denote by ti and ri respectively. Let Pi{t) be 
the probability that a particle occupies site i at time t: the master equation is then 

d 

— = ti+iPi+iit) ri-ipi-i{t) - (4 + ri)pi(t). (1) 
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For concreteness, we consider a periodic chain of N sites, but we are primarily 
concerned with propagation of particles on infinite chains: that is, we consider the 
limit of large N before any limit of large time. In this limit, propagators will be 
independent of the choice of boundary conditions. We also discuss finite chains with 
reflecting and absorbing boundaries in section [2TT] below. 

We use an operator notation where the ket |i) represents the state with the particle 
on site i, normalised so that = 5ij. Then, defining the state \P{t)) — '}2,iPii.'^)\^)^ 
the master equation can be written as ■^\P{t)) — VK(^^|P(i)) with 

N 

W^(i) = ^ - 1) + r,|z + 1} - {k + r,)K)] (*| (2) 
1=1 

N 

= ^(|z + l}-|j))(r,(i|-^,+i(» + l|). (3) 
1=1 

In this section, and wherever we consider systems with periodic boundaries, site 
i = TV + 1 is equivalent to site 1 and site is equivalent to site N . This allows 
terms to be rearranged as e.g. in going from ([2]) to ([31) above. 

We initially focus on an important special case: we assume that all rates are 
finite and that Wi^-i = Ili'"*- This ensures that all currents vanish in the long- 
time limit of the system. Under this assumption, we associate an energy with 
each site, measured downwards from an arbitrary baseline and determined through 
^i+ie^'+i — rie^'. Then, the model respects detailed balance with respect to the 
distribution p^*^ — e^' /{J^r^^"^)' where the sum runs over all sites. With this sign 
convention, a site with large positive Ei has a large Gibbs weight. The reason for this 
choice will become clear below. 

We now discuss general duality relations between pairs of hopping models. To 
this end, we re-parameterise the transition rates ii and by defining transition state 
energies E^_^_i associated with the links between sites: 

i.^er^'-i^^'\ ..=e-^^-^^^'^ (4) 

Note that site energies Ei and transition state energies E'j.). i are measured with respect 
to separate arbitrary baselines, so they do not necessarily have to be positive; but an 
interpretation in terms of activated hopping processes becomes problematic unless 
E, + Ei^i > 0. 

Two subclasses of these general hopping models have been studied quite 
extensively in the past [l[6l[i[10l[ll[T5l[16l[17]. The first subclass is the (pure) 
trap model, which is the case Ei_^_i = for all i, with the Ei being independent 
and identically distributed (i.i.d.); the second is the (pure) barrier model, which has 
Ei ~ for all i, with the Ei^i being i.i.d. A duality between pure trap and pure 
barrier models was discussed in [lOj . which we now generalise. The duality relation 
involves an inversion of the potential, which swaps the meaning of transition state 
energies and site energies: the various definitions are illustrated in Fig. [1] To present 
this relation, we define an alternative representation of the general hopping model, in 
which sites have indices i + 5 (throughout this article, i is always an integer). The 
hopping rate from i ~ ^ to i + ^ is Ri and the rate for the reverse process is Li. The 
master equation for this process has a similar representation in terms of state vectors 
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Figure 1. (Top) Illustration of an 'energy landscape', defined by supplementing 
the site energies Ei by transition state energies E-.i . The rate for hopping from 

2 

a site is related to the sum of the site energy (measured downwards as shown) and 
the adjacent transition state energy (measured upwards). (Bottom) On inversion 
of the potential, the transition state energies become energies of new sites with 
half-integer indices. Transition rates between these sites are controlled by the 
energies Ei which now have an interpretation as transition state energies. 



|i + i), with a master operator 

N 

W^(V2) ^ ^ (|, + 1) _ |, _ 1)) - i| - L,{^ + i|) . (5) 



i=l 



If we choose Ri = ri and Li = £i+i, then W^^^^'' describes the same model as 
W^^\ up to a simple relabeUing of sites (site i in W^^^ is mapped to site i — ^ in 
]/j/(i/2)^^ To describe instead the inverted potential of Fig. \T\ we choose 

Ri=e ■ 2 Li^c '+2 (6) 

With these definitions, the duality relation takes a simple form: we define 
operators 

N 

5= ^(K + l)-|z))(* + i|e-''-*-^ (7) 
1=1 

N 

= 5:|z)(e-^'-i(z-i|-e-^-i(^ + i|) (8) 



and 



N 



^ - E(N + 5)-N-5))(*|e"''' (9) 
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N 

= 5]|» + i)((z|c-^'-(z + l|e-^'+0- (10) 

i=l 

Hence, by combining the first expression ([7]) for 5 with the second expression pI7|) for 
5, and vice versa, one verifies that 

W^^^^SS, W^^'^^ =^S. (11) 

It foUows that 

51^(1) = (12) 

and that 

W'^^'^ S = SW^^'^\ (13) 

These two relations express the key duahty between the two master operators from 
which all our other results are derived. They imply, for example, that W*^^^ and M/^^/^) 
have the same spectrum of eigenvalues: if is a right eigenvector of PF^^/^) then 
S\ip) is a right eigenvector of W'^'^^ with the same eigenvalue. The exception is the 
singular case where S\^) — 0, which can happen only when the eigenvalue is zero. 
An analogous argument can be made in the other direction, showing overall that the 
nonzero eigenvalues and associated eigenvectors of VF*-^-* and VF*-^/^' are in one-to-one 
correspondence with each other. 

In fact, the structure of (jll[) occurs in the context of supersymmetric field theories 
(for a discussion of supersymmetry in the language of operators, applied to energy 
landscapes, see [18]; for a more general introduction, see [E]). These hopping models 
are very simple examples of supersymmetric partners: the spaces and {|i -I- ^)} 
can be interpreted as zero- and one-fermion subspaces of a generalised superspace. 
Acting to the right, the operator S annihilates all elements of and S annihilates 
all elements of {|i + i)} so S*^ = 5* =0: these operators are the supercharges of the 
theory. We can then write W^^^ + l^(i/2) ^ SS ^ SS which allows us to identify the 
two models as supersymmetric partners. 

It is also instructive to consider the continuum limit of our lattice model. Assume 
that the lattice spacing is a, so that the position of site i is Xi = ia. We define 
continuous functions Vi{x) and Vi/2{x) and take Ei — Vi{xi) and E^^^i — yi/2(2;i_|_i). 
Within the continuum limit, we represent the superspace explicitly, using a basis |x, n), 
where x is the position and n = 0, 1 distinguishes the zero- and one-fermion subspaces. 
That is, taking fermionic operators c, with = (c^)^ = and cc^ +c^c = 1, we take 
c|x,0) — and |x, 1) — c^|x,0). To make contact with the basis used for the lattice 
model, we identify |i) with \xi,0) and \i + i) with \xi, 1). If we then divide S and S 
by a and take the lattice spacing to zero, we arrive at 

5 = c-^e-^v.(-), ^^ct-^e-v-iW (14) 
dx dx 

The master (or Fokker-Planck) operators follow as 

l¥(i' ^SS^ cct Ae-^i/^(-)-le-^^(-), 
dx dx 

1^(1/2) ^55^ ctc-^e-Vi(.)_lg-v,/.(.) (15) 
dx dx 

This makes it obvious that the duality just swaps the trap and barrier parts of the 
potential, i.e. inverts the energy landscape. In the case without thermal activation. 
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which corresponds to V1/2 — ~Vi — V, the duahty reduces to the standard one for 
diffusion in the potentials V and —V. One part of this duahty was used recently in [50] 
to map boundary-driven steady states with current to current-free equilibrium states. 
Briefly, if P{x) is a steady state probability distribution of M^*^^-* in the presence 
of a boundary field so that W'^'^^P{x),0) = 0, then |P(x), 1) = S\P{x),0) obeys 
S\P{x), 1) = 0. Thus, \P{x), 1) is a steady state of W^^^'^\ but the current in this 
state is —e^ {d/dx)e~^ \P{x),l) = —e^c^S\P{x),l) — 0. The same approach also 
works on the lattice, and indeed the rate transformation (182) in Ref. [20] is the same 
as our 

The total master operator combining the dynamics in the zero and one-fermion 
spaces can be made Hermitian by a standard similarity transformation: with X = 
ccte5^i(^) _l_ ctce5^i/2(^) one has 

= cc^e~5^i(^) Ae-^i/2(^) Ae-5^i(^) + ctce"5^i/^("=)— e^^^^^)— e^^^i/^^"^) (16) 
dx dx dx dx 

Using cc^ = 1 — c^c to extract the extra contribution from the onc-fcrmion subspace, 
one can simplify this to 

dx dx 

^ (17) 



hvlixr - Vl/^ixf) kvl'ix) ~ Vl'/,{x)) 



For the 'standard' case with uniform diffusion constant V1/2 = —Vi — V, the one- 
fermion term simplifies to —c^cV"{x) and this is exactly the term that was used e.g. 
in [18[ I19j to construct dynamics that (in one dimension, and in the one-fermion 
subspace) converges to maxima rather than minima of the potential. 



2.1. Choice of boundary conditions 

In the previous section, we discussed systems with detailed balance and periodic 
boundaries. We now discuss how the duality relation applies on finite chains with 
reflecting or absorbing boundaries. We note in passing that some of these relations 
may be generalised both to boundary-driven models and to those with a finite bias 
acting in the bulk |2lj. However, for this work, we restrict ourselves to models without 
a steady-state current. 

In the case of periodic boundary conditions, the derivation of our duality relations 
required that the rates satisfy the global constraint Yli = Ili ; ™ order to guarantee 
detailed balance. An alternative that avoids this constraint is to use finite chains 
with reflecting or absorbing boundary conditions. These can be obtained by allowing 
zero rates in the periodic system. For a system with reflecting boundaries, we take 
(■I — riq — Q va. W'''^\ In the notation of energies, we set formally e ^ = 0: sites 1 
and N are then reflecting boundaries because they are separated by an infinite barrier. 
This has no effect on the equilibrium steady state, which now satisfies detailed balance 
whatever our choice for the remaining nonzero rates. The duality transformation 
carries through as before, resulting in a system with La? = i?i = 0. In this system 
there are no transitions out of site \ , which is therefore absorbing. The duality thus 
relates models with reflecting and absorbing boundaries, and as before these models 
share the same eigenvalues and their eigenvectors (for nonzero eigenvalues) are related 
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through the operators S and S. Note that for w'^^^'^\ the steady state that is reached 
in the long-time Hmit has the particle fully localised on the absorbing site \ . Detailed 
balance still holds, with e.g. the propagators obeying ([M]) below; the balance of 
transitions in the steady state is trivial because there are no transitions taking place 
at all. 

One can go further by distinguishing whether the absorbing site in 1^(^/2) is 
reached from the left (from site |) or from the right (from site N ~^), and accordingly 
split the absorbing site into two sites ^ and iV + ^ . The duality relation to the system 
W^^^ with two reflecting boundaries then holds as before. The only difference is that, 
because the system now has two absorbing sites, the master operator has two zero 
eigenvalues. The two corresponding right eigenvectors are localised on the absorbing 
sites, while the elements of the left eigenvectors give the probabilities that a particle 
initially on a given site will end up on either one of the absorbing sites (this will be 
discussed further in later sections). 

Finally, one may also take an operator W'^^^ with one reflecting boundary and 

-El 

an absorbing site at the other boundary: this can be done by setting e ' = 
and e~^" — 0, so that £i — In ^ = 0. In the dual model of this system, 
i?i = Rn — = 0, so that the reflecting boundary at site 1 maps to an absorbing 
boundary at site ^ and vice versa for sites N and N — ^. 

2.2. Alternative formulation, and generalised duality 

So far, our results apply to models in which the steady state has zero current (that is, 
periodic chains with global detailed balance, and chains with reflecting or absorbing 
boundaries). However, there is a slightly modified duality relation that holds even for 
periodic chains without detailed balance, where the steady state has finite current. 
Instead of factorising W'^'^^ = SS as in (fTTj) . one may instead write 

= ^DJ (18) 

withD = and J = N+5)(''i(*h^i+i(«+l|)- [In the continuum 

limit, one has W^^'' = —-^J where J = — e~^i/2(^) Ag-^iC^) jg the probability current 
operator, so the master (Fokker-Planck) equation has the form dfP = — V(JP) as 
usual, where JP is the probability current, which is linear in P. 
With these definitions, one may verify that 

t^(i/2) ^ _(ji5)t (19) 

or equivalently that 

W'^^^D ^ ~DJD ^ D{W'-^/'^^)'<. (20) 

If detailed balance holds, the relation {i + ^\W^^/^'>\j + ^)e^'+^ = {j + + 
i)e '+3 may be combined with (fT9|) to recover (fT3|) . (One way is to define E = 

J2i K + 5>(« + jle'^'+i and to note that D = -SE, J = E^^'S, and by detailed 
balance W^^^/^) = £:(PF(i/2))t^-i = -EJDE-^ = SS.) 

We also note that the operator D may not be inverted, since Dj^i N + |) = 0. 
However, if one considers for example a periodic chain of N sites with an absorbing 
site 1 and a refiecting barrier i, so that ri = £i = r^r = 0, then one may write 
= -D'J and 1^(1/2) = -{JD'Y with D' = D + In that case, it may be 

verified that the inverse of D' does exist, and the duality transformation becomes 

{D'y^W'-^^ D' = (I^(i/2))t (21) 
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Figure 2. Plots of Green's functions that illustrate the duality relation between 
the propagators. The plots were generated using the effective dynamics scheme 
described in section |4] but II26I I applies to propagation under the effective 
dynamics as well as to propagation with real dynamics. (Left) Propagator Gm,n, 
describing motion from site n to site m in a model of mixed trap- barrier type with 
Hb = 3 and fiq^ = 0.5 (see Sec. IS.ll for definitions). The 'squares' represent 'trap 
regions' within which the particle is locally equilibrated. The structure within the 
trap depends on the final site energies. These energies are randomly distributed, 
so that sites with low energy —Ei give rise to dark 'stripes'. However, local 
equilibration within the trap means that the propagator depends only weakly on 
the initial site n for all the n within a single trap. (Right) Propagator G n.m in 
the dual model of that shown on the left: note that m is now the initial site and 
n the final site. The probability for the final position is localised on a few deep 
'trap' sites. Eq. I|26|l states that for given values of m and n, the gradient of the 
right plot in the vertical (m)-direction is equal to the gradient of the left plot in 
the horizontal (n)-direction. 



We note that the transformation operators D and D' do not depend on the 
disordered rates and £i. An invertible transformation that is independent of disorder 
was used by Schiitz and Mussawisade [llj to study a reaction-diffusion model in a 
similar disordered environment. The duality transformation (enantiodromy) 2? used 
in [TT] is related to the transformation D' , and can be used to prove some of the 
relations for propagators that we discuss in the next section. However, analysis of 
the reaction-diffusion model is much more complex than the single particle system 
considered here. In the following, we restrict our analysis to the single-particle case, 
explaining which results may be obtained by alternative means. 

3. Propagators 

We now consider the propagators of these hopping models. For convenience, we first 
consider periodic boundaries in the case where detailed balance holds. Generalisations 
to reflecting/ absorbing boundaries are straightforward using the approach discussed 
above (Sec. 12. ip . For the model W^^\ we define the propagator Gn)m{t) as the 
probability that the particle is on site n, given that it was on site m a time t earlier: 

GilW = ("le'^''' V) (22) 
Similarly, for the model W^f^/^), we define 

G!^f,™+.W = (-+^N"'"*I-+^) (23) 
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Detailed balance holds for both models so we have, bearing in mind our sign convention 
for energies, Gnln{t)e^-^ = Gm^„(t)e^" or 

e-^-GW^W^e-^-GWjt), (24) 

with a similar relation for G'^^/^^. We now appeal to the duality relation SW'^^^ = 
to obtain {n + ij^Se'^'^V) = + \\e^^^'^''^'S\m) . Using the explicit form 
of S", we arrive at 

e--"G«„(i) -e--"-G,«,,„(0 = \c^g^^,At) - G^^t^-^^^)] ^^S) 

from which detailed balance for W'^^'^ implies 

GW„(t) - G«„^,(0 = G^f (t) - G^^f (0. (26) 

This result generalises Eq. (5) of Ref. [10]. Equation ([26]) may also be proved directly 
using the relation ([20|) . Hence, this relation between propagators applies on periodic 
chains, regardless of whether detailed balance holds. It holds for all times, including 
cases where the limit of large time is taken at fixed system size. 

Two comments are in order here. Firstly, a relation similar to (|26p was proven for 
propagators in the reaction-diffusion model of [11] . That model has pair annihilation, 
so we consider an initial state with two well-separated particles. If one then takes the 
limit iV — > oo at fixed time, one arrives at two independently propagating particles 
which can be shown to satisfy However, the mapping of [TI] cannot be used to 

prove for large times and periodic chains of fixed length: in that case, it would be 
natural to use an initial state with exactly one particle but the mapping then breaks 
down (except for specific cases where the periodic chain is broken into segments by 
absorbing sites and reflecting barriers). Secondly, for the specific case where the Ei 
are chosen freely but E^j^^ = for all i, the relation ([26]) was given in Ref. [10], but it 

is also implicit in Ref. [22]. (In the notation there, it reads f[q)^^q'{z) — T'^^, (z) f {q' ) 
and can be derived by inserting Eqs. (2.15) and (2.16) into (2.19) of Ref. [22] and 
expanding appropriately.) 

The equality relates differences with respect to the initial site of propagators 
in the two hopping models. If the propagator for one model is known, the other 
can be calculated by successive application of this equation. An illustration of this 
relation is shown in Fig. [2] We also note that ([26| is symmetric between the two 
master operators: any set of transition rates can be interpreted either as a model of 
type W^"^^ or as a model of type VF^^/^). To reinforce this symmetry, we note that the 
relation (n|S'e^''^"*|m - ^) = (n|e^*''*5'|m - i) implies that 

^^^-iG^/^^^_^(t) - e^'^-iGi^^)^., (t) = [Gl%{t) - Gi^LiW] e^^'-i (27) 

which is the same relation as (|25p . but with the original model W^^'' expressed in the 
form and vice versa. 

The simple form of pS)) means that we can take the disorder average. As 
long as the disorder is translationally invariant, e.g. if the transition state energies 
{Ej^_^_l} and site energies {Ei} are taken from two (possibly different) translationally 

invariant distributions, the disorder-averaged propagators G„.m(i) will depend only 

on k = n — m. Defining then Afe = G'l^\t) — G^^l^\t)^ the relation becomes 

A_fe = A_fe_i. Thus, Afc is a constant; but G''j^\t) = 1 and similarly for G']^^'^\t), 
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so A^; = and the constant has to vanish. One concludes that G[^''(t) — G^^^j^\t). 
Re-instating the notation with separate initial and finite site labels, we have 

G'L!n(0 = g\^Li [i) = G^mji {i) (28) 

where the last equality follows from left-right symmetry. Remarkably, then, the 
disorder-averaged propagators of the dual models are equal on all scales of length and 
time. Again, this generalises our earlier statement |10| relating average propagators 
in pure trap and barrier models. A similar result applies for disorder distributions 
which break left-right symmetry, as long as they remain translationally invariant. 
Take for example a model described by W^^^ with site and barrier energies from a 
translationally invariant distribution, and then modify the rates according to a site- 
independent prescription, so that the left-going and right-going rates acquire different 
distributions. With periodic boundary conditions, translation invariance again holds 
and we obtain by the same arguments as above 

GS„W = Gi^™ W, (29) 

with both sides again depending only on the difference k = n — ni. The same result 
applies on long chains with reflecting boundaries, as long as we are far away from 
these boundaries so that translation invariance is not broken. 

4. Effective dynamics scheme 

Several effective dynamics and renormalisation schemes have been proposed to 
approximate propagators for motion in random potentials [HI [T31 [T71 [53] . The most 
notable success in this area is the work of Ic Doussal, Monthus and Fisher [T^] , which 
we refer to as DMF. They considered the Sinai model [H], in which the and ii are 
independently and identically distributed, so that the site energies Ei follow a random 
walk in real space. For long time scales in that model, DMF found a renormalisation 
group (RG) scheme that can be implemented analytically and gives exact results for 
a variety of observables. 

Later, Monthus ^17j applied a related scheme to the pure trap model, in which 
(■i = Ti and these rates are independently and identically distributed with a power- 
law distribution P{ri) oc r'f'^^'^ ^. We refer to this scheme as Mon03: it can be 
treated analytically and its predictions are exact in the limit of large /i (our notation 
follows that of [TUj: to arrive at the notation of [T7], replace fi by i). More recently, 
Monthus and Garel (MG) derived a general RG scheme [53j that is not restricted to 
one dimension nor to single-particle systems. Like the schemes of Refs. [TH[T7], the 
MG scheme is effective when disorder is the dominant source of fluctuations. Finally, 
we recently |10| introduced a modified scheme for pure trap and barrier models that 
respects the duality symmetry ([26]) but requires a computational implementation. 
This prescription accounts for effects that were neglected in previous schemes, and 
reduces to the Mon03 scheme in the limits in which that method is exact. 

Here, we generalise the scheme of [10] to general disordered potentials (in 
one dimension). We arrive at a method that encompasses the DMF and Mon03 
methods in the limits when they are exact, and is also consistent with the duality 
relation ([26]) . Our method also shares features with that of MG, but it respects the 
duality symmetries discussed above, while the MG scheme does not. We describe the 
application of this scheme to models parameterised in the form of W'^^\ with site 
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energies Ei and transition state energies ■ For notational convenience, we simply 
denote these master operators by W from now on. 

Definition of the effective dynamics scheme 

The idea of the effective dynamics scheme is to describe motion on long time scales in 
terms of a coarse-grained set of co-ordinates in which the system evolves slowly. We 
have developed two versions of the scheme: in the main text we describe the more 
physical and intuitive version that also gives the most accurate description of motion in 
the system. We discuss the justifications of our scheme in later sections ( [Appendix A| 
and Sec. [5]): in particular, the second version of the scheme can be justified more 
formally, and it gives similar results, but with larger errors for finite 

The scheme is parameterised by a time scale . On that time scale, we partition 
the Id chain into effective trap and effective barrier regions. The coarse-grained slow 
co-ordinates are the occupancies of the effective traps: 



Po.it) = £ P^{t)■ (30) 



We continue to use Roman indices i,j,... for the original site labels of the model, 
while Greek indices a, (3, . . . label the slow co-ordinates. Thus, Oq and are the 
leftmost and rightmost sites within effective trap a. The traps form an ordered set 
along the chain, so we have ba-i < Oa < ba < Oc+i. We refer to the regions between 
these effective traps as effective barriers. 

For a given value of F, the effective dynamics scheme gives (i) an approximation 
for the propagator Gmn(r^^) and (ii) an approximate equation of motion for the slow 
co-ordinates, valid for times t > T~^. In addition, the scheme specifies (iii) how the 
set of effective traps evolves as F is decreased towards zero. 

Physically, the idea is that motion within effective traps is fast, while motion 
between traps is slow. Assuming equilibration within effective traps, we have for sites 
within an effective trap region, 

G,,(F-i) ~ e^^-^°, a„ < i,j < 6„ (31) 

where we identify the free energy of effective trap a: 



e^" = Y: e"^- (32) 



If site i is between traps a — 1 and a, that is &q-i < i < Oq,, then it may relax into 
either trap before equilibrating there, and we have 

' (l-t;,<"~^V^-^°-S aa-i<j<&«-i 



aa < j < ba 



G,.(F-)^ ; ' ""-^"^f (33) 

[ vl ^ e^3 

where wj" is the probability that a particle initially on a site i between traps 
a — 1 and a relaxes first into trap a. These probabilities can be obtained by 
considering propagation on a chain with absorbing sites and Oq, and calculating 
the probability of absorption into each of these sites for a given initial site i. The 
result is 

^ = e 2 > e '+2 (34) 
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with 

e^-h = e^'^+h. (35) 

This completes point (i) above. Turning to point (ii), the equations of motion for the 
slow degrees of freedom are 

dtPa{t) = ta + lPa+l{t) + ra-lPa-l{t) - + ra)Pa{t) (36) 

with 

Tq = e -^Q = e (37) 

That is, hopping takes place only between nearest neighbours, and the rates are given 
by the free energies of the effective traps and 'transition state free energies' F^_^_l. 

It remains to discuss point (iii) above: how the set of effective traps depends 
on the time scale F. Initially each effective trap contains a single site i, and each 
effective barrier a single transition site i + \- Each stage of our effective dynamics 
begins by calculating rates for motion from trap a: we calculate a renormalised rate 
pa associated with the rate r^, i.e. with leaving the trap to the right: 

Pa = + +2r„(^„+£„+l) + (C-^„+l)2. (38) 

and similarly for motion to the left 

Aa = ^(^a + r-a + r„_i) + 1^^2 + 2ia{ra + Ta-i) + {r^ ~ r„_i)2, (39) 

Notice that pa is symmetric in and la+i- In a (renormalised) landscape similar to 
that for pure traps, would be comparable to r^ and so pa ~ if the other rate 

{(-a+i) is small. Conversely, in a landscape resembling pure barriers, one would have 
£q,+i of the same order as and so to leading order pa ~ r^ + ia+i- The expressions 
(|38l39p cover both of these limits but extend them to general landscapes. 

Having calculated the pa and Aq,, we select the largest of these rates across all 
traps, and update F to this largest rate. At this point, several cases arise, which are 
illustrated in the following section, with reference to the example landscape in Fig. [31 
Here we give the rules: Supposing that the largest rate is pa-, we now remove trap a. 
We make a case distinction, depending on whether the local landscape is nearer to the 
pure trap or the pure barrier case, as in the discussion of the rate-dependence of pa 
above. If la+i > (-a we combine trap a into trap a + 1, leading to a new effective trap 
containing sites Oq . . . 5q.+i. Barrier region ^ is removed. Conversely, if (-a+i < ^a, 
we remove trap a from the list of slow co-ordinates: this amounts to combining barrier 
regions a ± i . For the case where the largest renormalised rate is Aq , associated with 
motion to the left, the rules are similar, in accordance with left-right symmetry: if 
Tq-i > ra we combine traps a — 1 and a; otherwise we remove trap a which amounts 
to combining the two barrier regions a ± i. 

Having combined the appropriate traps or barriers, we finally recalculate the free 
energies Fa and F^^i that are affected by the change, and hence obtain new hopping 
rates ra and ia, and new renormalised rates pa and Aq. From here on the process is 
iterated, i.e. we find the largest rate among the {pa,^a}, update F, and merge the 
appropriate traps or barriers. 

For any given F, we can therefore calculate the approximate propagator and the 
equation of motion for the slow degrees of freedom. We emphasise that while we have 
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Figure 3. Illustration of the RG scheme introduced in this paper, discussed in 

the main text and in [Appendix A| As in Fig. [l] sites are shown with closed 

circles and transition states with open circles. However, we now introduce 

effective trap regions, represented by filled ovals, and effective barrier regions, 

represented by open ovals, (a) Energy landscape for a hopping model with rates 

satisfying l|40|l . The fastest rate r\ is associated with the energy difference 

E\ + E3 : the link associated with this rate is emphasised by a bold line. 
2 

(b) Energy landscape obtained after one stage of renormalisation, with slow 
degrees of freedom {. . . ,po(t)i P2{t), ■ • .}• The transition states ^ and ^ and 
site 1 have been incorporated into a composite transition state, which we refer to 
as an effective barrier region. Within the effective dynamics, a particle initially 
on site 1 relaxes either to site or site 2, according to II33I I. The link associated 
with the largest rate in the renormalised master operator is again emphasised. 

(c) Energy landscape after the second stage of renormalisation, with slow degrees 
of freedom {. . . , po(i) + + P2(i), . . ■} The two effective traps on sites and 
2 have been combined into a single effective trap. Within the effective dynamics, 
particles initially on any of these sites relax to an equilibrium distribution over 
sites 0, 1, 2. 



defined rates ra, £a, Xa, Pa, free energies Fa,F^_^_i and parameters aa,ba,v^ ^ , all 
of these quantities are fixed if we specify the set of slow co-ordinates (parameterised 
in terms of the Oq and ba), together with the full set of bare energies Ei^E^j^i. We 
discuss below how the progress of the scheme can be written in terms of a F-dependent 
projection operator, using the notation of Seed] However, we first give an illustrative 
example of the effective dynamics in action. 

4-2. Example of effective dynamics and comparison with DMF scheme 

We illustrate the effective dynamics using the example landscape shown in Fig. [Sj 
concentrating on sites 0, 1, 2 of a long chain. Consistent with Fig. [5^, we choose the 
rates to lie in three well-separated sectors: 

ri,^i >ro,4 >4,r2. (40) 

For concreteness we also take ri > li (i.e. E\ > -E3), I2 > rn and E2 > En. 

2 2 

At the initial stage of the dynamics, sites 0, 1, 2 each constitute an effective trap: 
the effective equation of motion ([36]) coincides with the original master equation 
(fTj) and the relevant degrees of freedom are simply the original site occupancies 
{. . . , pQ{t) , pi(t) , p2{t) , . . .} where the (...) indicate that we are concentrating on part 
of a large chain. Taking the original rates and £i, we evaluate the parameters pa and 
Xa, and the largest of these is pi = ri+£i+e2 + ^\/ri + 2ri{£i + £2) + {ii - l-iY ■ (To 
leading order in the largest rates r\ and p\ = \\ — r\^ £\\ but the correction from 
£2 can then be shown to make p\ > Ai.) Thus, the first step of the effective dynamics 
occurs as F is decreased through pi. From (f40|) . we have £2 < £1, so the rules state 
that we remove the trap on site 1. Physically, the idea is that the transition state 
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energies Ei and Es are similar to each other, so they are combined into an effective 
barrier. On the other hand, the particle spends very little time on site 1, compared to 
sites and 2, so the co-ordinate pi{t) relaxes quickly to a small value and is no longer 
a relevant (slow) co-ordinate. 

The remaining slow co-ordinates are therefore {. . . ,po(i)jP2(i), ■ • •}, and the 
energy landscape on this time scale is shown in Fig. ^p. From (|36p . the equations 
of motion for these two co-ordinates at this stage are 

dtp^{t) = r_ip_i(t) +4p2(t) - {r'o+iMt) 

dtP2{t) = r'^P^it) + tsPsit) - {r2 + Qp2{t) (41) 

where the unprimed rates and are the original hopping rates among the sites 
of the model, but two new rates have appeared: Tq = e"^""^ and (.'2 = e~^^~^. 

Here e^ = e^i + e^i , and F is the 'transition state free energy' for the effective 
barrier region between sites and 2, constructed according to (|35p . One may also 
construct the effective propagator on the time scale in accordance with ([3T1 [33]) : 
the non-zero matrix elements among sites 0,1,2 are 

Goi~e^i"^, Gsi-e^i""^, Gqo = G22 = 1- (42) 

where the approximate equalities simply indicate that these are propagators under the 
effective dynamics. 

For the next stage of the effective dynamics, the largest hopping rates are rg and 
£'2- Because E2 > Eq, the former will be larger than the latter, and the remaining 
rates £0 and r2 are much smaller by our assumption ((40|) . Correspondingly, when 
we calculate rates p and A for the next stage of the effective dynamics, the largest 
one win be p'„ = ^{r'^ + £0 + £'2) + \^Jr2 + 2ro{£'2 + £0) + {£'2 - £oY ■ From the 
assumptions (|40|) . we have £'q <C £2, and the rules of the effective dynamics state 
that we merge sites and 2. Physically, the site energies Eq and E2 are similar so 
these are combined into an effective trap which also contains the intervening site 1. 
In accordance with (|30p. the slow degrees of freedom are now simply {. . . ,pm2{t), . . .} 
where poi2(0 = Po{i) ~^ Pi{t) + P2{t) is the occupancy of an effective trap containing 
sites 0, 1 and 2. The corresponding energy landscape is shown in Fig. [31:. The equation 
of motion for poi2(i) at this stage is 

dtP012{t) ^ r-ip-l(t) +£3P3{t) - (?-012 +4l2)P012(i) (43) 

where the new rates appearing at this stage are roi2 = e 3 for hops to the right 

— ^ ^ 012 — 'Y 

from the effective trap, and ^012 = e "2 for hops to the left. Here, F012 is the 

trap free energy, e''^"'^^ — + e^^ + e^^, in accordance with ([32]) . The approximate 
propagator on time scales l/p'o can be constructed from ([3T | [33 l) . giving 

Gj, = e^^--^oi^ , 0<i,j <2. (44) 

which is independent of the initial site i as long as it is within the trap, consistent 
with the idea of local equilibration. 

4-3. Comparison with other schemes 

It is useful to compare this scheme with the effective dynamics of DMF [T^]. In 
that scheme, the slow degrees of freedom Pa{t) are simply a subset of the original 
site occupancies Pi{t), so each effective trap contains exactly one site: Oa = ba for 
all a. Further, instead of calculating transition state free energies as in ([SS]) . one 
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Figure 4. Illustration of the RG scheme of DMF, discussed in the main text, (a) 

Energy landscape associated with a master operator satisfying II40I I. The fastest 

rate ri is associated with the energy difference Ex + E3 : the link associated with 

2 

this rate is emphasised in bold, (b) Energy landscape associated with the model 
after one stage of renormalisation using the effective dynamics of DMF. The fast 
motion on these time scales involves a particle initially on site 1 relaxing onto 
site 2, with probability unity. The link associated with the largest rate of the 
resulting master operator is again emphasised, (c) Energy landscape after two 
stages of renormalisation following DMF. Particles initially on sites 0, 1 or 2 will 
relax onto site 2. 

takes simply e^°"i — maxfc^_j<j<Q^ e ^+^, assuming that the transition state free 
energy of each barrier region is dominated by the largest transition state energy within 
that region. Finally, one also takes simply Ac = (-a and pa = ra when deciding 
which traps to remove as F is reduced. The progress of the scheme is illustrated in 
Fig. m for the same example landscape considered in Fig. [31 The first step of the 
renormalisation scheme takes place when F = ri and the slow degrees of freedom 
after this step are {. . . , po{t) , p2{t) , . . .} as in our scheme. In the DMF scheme, the 
equation of motion for these degrees of freedom is of the same form as (|4ip . but the 

new rates are r'g — e and i'2 — g i [note that Ei > Es, from ([3D]), which 

sets the transition state free energy]. Further, in the DMF scheme, all elements of 
the approximate propagator are zero or unity: for times of order 1/ri, the non-zero 
elements of the approximate for < j < 2 are 

Goo — G21 — G22 — 1- (45) 

which can be compared with (|42)) . 

Then, the second stage of the scheme takes place at F = rg, after which the slow 
degrees of arc {. . . ,p2(t), . . .}, and the new rates for motion from the trap containing 

site 2 are £'2=6 ^ and = e ^ ^ . For times of order I/tq, the non-zero 

elements of the approximate propagator between sites 0, 1, 2 are 

G20 = G21 = G22 = 1 (46) 

which can be compared with (|44p . 

In general then, the DMF scheme gives different results to our scheme, although 
broad features are similar. In the limit where all rates are well-separated: ri ^ £1 ^ 
I2 ^ To ^ £Q,r2, it may be verified that the two schemes coincide. That is the 
condition in which the DMF scheme is exact, and we conclude that our scheme is 
also exact in that limit. In fact, our scheme is also exact in the less- restricted limit 
of well-separated rates given in ([¥0)1 . while the DMF method is inaccurate in that 
case. The clearest differences between the schemes occur in the propagators. For 
example, compring with pS)) . the DMF scheme ignores the possibility that a 
particle originally on site 1 may relax onto site at this stage, while such transitions 
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happen with probabihty approaching i if ri k, Ix. AdditionaUy, under the same 
condition ri ~ the transition rates Tq and £3 differ between the schemes by a factor 
close to 2. Compared to DMF, these rates are smaher in our new scheme, reflecting 
the possibility that a particle that hops from site to site 1 may return to site before 
visiting site 2. As discussed in [10], factors such as these must be included in effective 
dynamics schemes in order to obtain the correct scaling behaviour for models where 
many of the E^j^i are approximately equal. 

We also compare the scheme given here with that of MG [23 . That scheme 
states the slow co-ordinates at each stage and their equations of motion, although 
the explicit propagator Gmn among the original sites of the model is not given. In 
fact, for onc-dimensional hopping models, the 'full' MG scheme reduces to a simplified 
version of our scheme, in which effective traps are always simply removed from the list 
of slow degrees of freedom, but traps are never merged. In the language of Sec. 14.11 
one always assumes that the landscape has a 'pure trap' character. In choosing which 
traps to remove, one takes the parameters \a = Pa = fa -\- (-a, consistent with that 
assumption. One may verify that this procedure is exact and coincides with our 
scheme in the limit where all site energies Ei are well-separated from each other: it 
is therefore appropriate for systems such as pure trap models. However, the scheme 
does not preserve the duality of Sec. [2| in particular, while it is appropriate for pure 
trap models, it fails for their duals, which are pure barrier models. 



^-4-. Operator representation of effective dynamics 

The effective dynamics scheme can be interpreted as a projection of the master 
operator W onto its slow degrees of freedom. Briefly, any master operator can 
be diagonahsed as = ~^\\^r)^{^l\ and its propagator written as e^* = 
X^A I^-r)c Time scales are (globally) well-separated if there is a time t such 

that all of the exponential factors are either negligibly small or close to unity. For 
times t with that property, the time evolution operator is well-approximated by a 
projection operator 

e^^^T'oxr^)^ (47) 

A<t-i 

which can be rearranged into the form e^* k, \Pa){(la\ where the states \pa) and 
{qa\ have non-negative elements [25]. For the specific case of hopping models, the 
states \pa) indicate the effective traps into which the system relaxes on the time scale 
t, while the vectors {qa\ indicate the probabilities of relaxing into these states. 

For the effective dynamics, the key point is that motion on time scales longer 
than t = can be well-described by a renornialised master operator WR.ox(r) = 
Pcx(r)W^7'cx(r), since eigenmodes with At ;2> 1 that are irrelevant at time t are also 
irrelevant for all longer times [more precisely, e^* = e'^'^ "'''-^''* -I- 0(e~^*) and the 
error becomes small for t ^ F"^]. In the effective dynamics a set of effective trap 
and barrier regions corresponds to an operator ViT) that approximates Pex(r). We 
therefore define the renormalised master operator 

Wk{T) = V(T)WV{T) (48) 

and we note that a corresponding approximation for the propagator is 

G,„„(r^i) ~ (m|e^«/r|n) ~ {m\V{T)\n) (49) 
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For consistency with (j36p above, the renormalised master operator should take 
the form 

W^R(r) = 5](|Pa+i) - |Pa))(r„(Q„| - tc.+i{Qc.+i\) (50) 

a 

where \Pa) and {Qa\ are vectors associated with trap a. Thus, the effective dynamics 
represent a renormaUsation scheme in the sense that the operator Wr maintains the 
same form as the original master operator ^ as the cutoff F is reduced. To obtain 
such a form for PFr , we take 

P(F) = ^|P„)(Q„| (51) 

a 

where the \Pa) and {Qa\ are approximations to slow eigenvectors of W . Our scheme 
as given by Equs. (|30ll37p corresponds to the choice 

{i\Po.) = e^'-^° (52) 

and 

V- , <i<aa 

1, aa<i<ba (53) 

l-V^ ^ , ba <l< Oa+l 

In [Appendix A[ we show that the \Pa) and {Qa\ constructed in this way are indeed 
good approximations to slow eigenvectors of VF, under conditions discussed below. 

Thus, the effective equation of motion (I36p corresponds in operator notation 
to the equation dt\P{t)) = WR{T)\P{t)) = 'PiT)WP{T)\P{t)) and the approximate 
propagator of (|31|33p corresponds to G'mn(F~^) {m\P{T)\n). The quality of these 
approximations depends on two considerations. Firstly, the validity of the effective 
equation of motion p6p depends on the extent to which the projection operator 'P(F) 
approximates VcxO^)- Then, the extent to which the propagator G'mn(r~^) coincides 
with (m|P(r)|n) depends in addition on a separation of time scales, as can be seen 
from the discussion of (|T7)l . 

We observe that if 7'(F) w 'Pcx(r) at all stages in the scheme then the errors 
associated with the effective dynamics remain small as F is reduced, while large errors 
arise if 'P(F) becomes different from 7-'cx(r)- The conditions under which the operator 
V{T) represents a good approximation to 7^cx(F) are discussed in [Appendix A[ Some 
numerical tests are also given in Sec. [S] We summarise here the analytic results 
of [Appendix A[ On long time scales (small F) we identify the 'fastest relevant rates' 
ia and Tq, which are comparable to F. For a consistent renormalisation flow, we require 
that these fast relevant rates are typically much larger than all other relevant rates in 
their neighbourhood, except that (i) fast relevant rates Tq, may be comparable either 
to ia or to ia+i, and (ii) fast relevant rates ia may be comparable either to r^ or 
to ra-i- (The condition of globally well-separated time scales described above is not 
required: it is sufficient that eigenvalues for motion in the same spatial neighbourhood 
should be well-separated.) 

4-5. Renormalisation and duality 

To conclude this section, we discuss duality relations for the effective dynamics, 
restoring superscripts to distinguish between VF'-^-' and W^^^^\ If we renormalise 
an operator W^^^ = SS according to our scheme, we arrive at a renormalised model 
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that can be written in the form = SykSyk,- The operators S'r and S'r have the 

same form as S and 5, except that sites i are replaced by effective traps a, transition 
states by effective barriers, and energies by free energies. An important property of 
our RG procedure is that if we apply it to the dual master operator Vt^f^/^) = SS, 
we find that this renormalises precisely to VF^^^^^ — S^S^. Thus, in addition to 
the basic requirement that M^R(r) takes the same form as W, our scheme also obeys 
the general duality relation under landscape inversion. This is of course desirable, as 
renormalisation schemes should respect all symmetries of the models of interest. 

In the illustrations of Figs. [3] and HI the duality property follows for both the 
new scheme and that of DMF, because acting on these illustrations with the inversion 
operation of Fig. [T] leads to the same renormalisation flows that would be obtained by 
starting with the dual of the original model. 

Mathematically, the duality can be shown as follows. At each stage of the RG 
flow we have effective trap regions associated with wj^\ with intervening barrier 
regions. Assigning integer indices a to the traps, the barriers can be associated with 
indices a + i. In wj^^^^ the regions with integer indices a become effective barriers, 
while those with indices a + ^ become effective traps. The associated free energies 

Fa and F^+^ are the same in both cases. In the dual model W^^'^\ the rates for 
hopping to right and left from trap en — ^ are r^_i = ia and i^-^ = ^a-i, consistent 
with One also easily checks that the rates A^_|_i and Pa+^ that are obtained on 
renormalising W^^^'^'' are the same as those obtained on renormalising W''^\ according 
to Pct-i = Aq and Aq_|_^ = pa- Thus, supposing that we combine traps a and a + 1 

in an RG step on the model we also remove barrier a + ^. In the dual model 

we remove the trap with index a + ^, which corresponds to combining the 
barrier regions a and a + 1. Finally, it can be verified that be verified that the rules 
for merging and removing traps do preserve the duality between wj^'^ and wj^^^\ 
A brief comment is in order on the construction of S'r and Sr. From wj^'^ = 
and W^'^^ = SS one might naively identify Sr = V^^^S, Sr = SV^^^; 
but this choice does not satisfy the duality requirement that W^^^'' = SrSr. A little 
thought shows that one requires instead Sr = pWsV^'^/'^'> and = T'<^^/'^^S'P^'^'>. 
Here P^^/-^) = J2a \Pa+^){Qa+^ \ is the projector onto the effective traps of the dual 
model, with |-Pq+ i ) and {Qa+i I constructed in the obvious manner using the duals of 
(|52p and (|53|) . The operators Sr and Sr defined in this way are indeed the effective 
trap and barrier analogues of ((7 HT0|) . and the desired dual expressions = SrSr 

and wj^^^'^ = SrSr therefore hold. We note finally that not only the master operators 
but also the propagators produced by our RG scheme obey the required duality: the 
propagators are the matrix elements of the projectors P*^^^ and 7^(1/2)^ and one verifies 
by direct calculation that (m|p(i)(|n) - |n+ 1)) = {n + ^\V'^^/'^'> {\m + ^) - |to - i)) 
in accordance with the general duality relation (pSj). 

5. Specific disorder distributions 

5.1. Mixed trap-barrier models 

The renormalisation scheme that we have discussed can be implemented 
computationally without undue difficulty. It allows rapid estimation of propagators in 
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these hopping models, both for fixed disorder and for disorder-averaged properties. We 
first consider a model obtained by mixing the pure trap and barrier models defined 
above. In the pure trap model, the transition state energies are i^j+i — 0, while 
site energies are chosen from an exponential distribution with a mean of /x; that is, 
P{Ei) — {1/ fi)e~-^'/^^ with Ei > 0. In terms of rates, this implies ii — ri = Wi 
with P(wi) = {1/ ^)wl^^^^ ^ for < Wi < 1. Similarly, pure barrier models have 
Ei = for all sites, and transition state energies are exponentially distributed: that 
is, Li — Ri = Wi with the same distribution P{wi). 

We mix these models by taking both site and transition state energies to be 
exponentially distributed with means /ix and /iB respectively. The dynamical scaling 
of these models therefore depends on the parameters fiT and /iB • For the pure barrier 
model (^T 0), sites have half- integer indices and moving a distance r typically 

E 1 _ 

requires the crossing of a barrier i whose hopping rate is = e '+ 2 ^ r . The rate 
for actually crossing this barrier is suppressed because the particle is delocalised in an 
effective trap whose width is of order r. In the language of the effective dynamics, the 
landscape consists of wide effective traps separated by isolated transition states, and 
each site i within the trap contributes contributes e'^' = 1 to e^" . The result is that 
the time taken to move a distance r is T(r) ~ e ° ^ r/wi ^ r^^^^. In the pure 

trap model, the typical time for escaping from sites with large Ei is 1/wi ~ e^% but 
the barrier regions on this time scale are typically of width r and their free energies 

F 1 

therefore also scale as e ~ r, reflecting the probability of reabsorption in the 
original trap before arriving at a new one |10j . Thus, the typical relevant time scale 
is again T(r) ~ r^+''. 

We define the dynamical exponent z through the relation t ^ r^' , and identify 

z = 1 + /iB, /^B > 1, Mt ^ (54) 

with an analogous relation if /j,t > 1 and /iB ^ 0. 

In the mixed model, the crucial case distinctions are then whether the /^b and 
/XT are larger or smaller than unity. For example, if /^t < 1, the average value of 
e^^ is finite, and the free energies of relevant trap regions scale as e''^' ~ r. On the 
other hand, if /ix > 1, the site-averaged e'^' is no longer finite, and the sum in ([32|) 
is dominated by the largest site within the effective trap. In this case e''^' ^ r^"^ . A 
similar argument applies for /iB- Combining these results, we arrive at the dynamic 
exponent for the mixed model 

z = max(l, /iB) H- max(l, /^t) (55) 

which reduces to the pure trap case if /iB = and the pure barrier case if /ix — 0. Fig. [5] 
shows numerical results that are consistent with ([55)1 . Thus, the effective dynamics 
provide a natural framework in which to derive this kind of scaling result, although the 
above predictions for the dynamical exponent could presumably be obtained by other 
means. We also note that the renormalisation arguments of DMF give z — hb + fJ^T 
which is the correct result when both fiB and /ix are greater than unity. This is 
consistent with our assertion above that if fj,B , fJ-T > 1 the free energies of effective 
traps and barriers are typically dominated by single sites, and this is precisely the 
limit in which the scheme of DMF is valid without approximation. To be precise, it 
follows from the discussion of [Appendix A| that our effective dynamics scheme is exact 
in the limit where either /^b or /^t is very large, while the DMF scheme is exact in the 
limit where both /ib and /ix are very large. 
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Figure 5. (a) Mean square displacement, sliowing real dynamics (symbols), 
effective dynamics (solid lines) and power law predictions of l|55|l (dashed lines). 
The DMF scheme (not shown) also gives the correct scaling as long as > 1 

but gives the wrong exponent for the case where /^x = 0.5. (b,c) Propagators for 
real dynamics (b) and effective dynamics (c) for the case /^b = 10 and /.ix = 1- The 
time is t = 2^^: the propagator for the real dynamics is obtained by simulating 
its dual model and using l|26| l. since models with large /^x can be simulated 
much more efSciently than those with large ^b- Since /^b is large, nearby rates 
are well-separated and the effective dynamics gives a good approximation to 
the propagator even for fixed disorder. However, there are deviations in some 
neighbourhoods which are associated with the presence of eigenvalues of the same 
order as 1/t, as discussed in the text. (d,e) Propagators for the formal effective 
dynamics scheme (see [Appendix A\ and the DMF scheme (see Sec. 14.21 1. The 
formal scheme differs slightly from the physical scheme: for example, the formal 
scheme gives dj = for i = 131, 132, independent of j, while the physical scheme 
gives small finite values that are more consistent with the real dynamics. In the 
DMF scheme, all effective traps consist of single sites, so the scheme does not 
capture the broad 'effective traps' that are visible in the real dynamics. 
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In Fig. IH^bjC), we show example propagators obtained using real and effective 
dynamics for the case /ie — 10, /ix — 1- For these parameters, time scales in a 
given neighbourhood are sufficiently well-separated that the effective dynamics gives 
a good approximation to the propagator. In this case, it appears that the largest 
deviations between real and effective dynamics come from neighbourhoods in which 
there is an eigenvalue of W of the order of \/t. In the discussion of Sec. 14. 4| we noted 
that such deviations are expected even when our scheme gives P{T) exactly equal 
to Pox(r), and that the specific deviations seen at any time t should decay as time 
increases, so that the exact and approximate propagators remain close. Our results are 
consistent with this expectation. For example, comparing Figs.[5jD and[5j;, the effective 
dynamics indicates that sites in the vicinity of site 110 are separated into three traps, 
containing sites 106-113, 114-118, and 119-121: the effective barrier regions are simply 
single transition states. However, the real dynamics reveals that initial sites 114 — 118 
propagate both within that effective trap, and into the adjacent traps. On time scales 
much shorter than one would expect localisation within this trap; on longer time 
scales the trap will either merge with an adjacent trap, or become incorporated into 
an effective barrier region. We conclude that we have measured the propagator during 
this crossover, the details of which are not captured by the effective dynamics. 

In Fig. [Hdje), we compare the physical effective dynamics of l4.1l with the formal 
scheme of [Appendix A| and the scheme of DMF. As discussed in Sec. 14.31 the DMF 
scheme assumes that all effective traps consist of only a single site, so for any initial 
site i, there is a single final site j such that Gji — 1, with Gji = for all other 
final sites. It can be seen that the DMF scheme does identify final sites j with large 
Gji, but it underestimates Gji for other j. As discussed in [Appendix A[ our physical 
effective dynamics scheme means that for a given «, Gji is finite for j within contiguous 
regions of the chain; on the other hand, the formal scheme leads to Gji that is finite 
on a restricted set of sites within such regions. As discussed in the caption to Fig. [5l 
the formal scheme therefore underestimates Gji for some sites j. Since /ie = 10 
is quite large, we expect both of our effective dynamics schemes to mimic the real 
dynamics quite accurately, consistent with the data. Deviations between the two 
schemes and the real dynamics would increase as is reduced (data not shown). 
Similarly, the predictions of the DMF would mime the real dynamics more closely if 
/zt were increased since that scheme requires both /ix and /xb to be large. 

Moving to disorder-averaged properties, numerical results indicate that the long- 
time behaviour in these systems is associated with a scaling form of the diffusion front, 
as expected. That is (with k the distance between initial and final site as before), 

Gu{t)~<J-^g[k/a). (56) 

where the function g{x) is independent of the time t and we use a — a[t) — l/Go{t) as 
an estimate of the length scale associated with motion on a time scale t, ensuring that 
g{0) = 1. Our numerical results then indicate that the shape of the diffusion front 
g{x) depends quite strongly on max(/iB, /ix) and much more weakly on min(/iB, /ix)- 
As in Ref. [10], the effective dynamics give good agreement with the real dynamics 
when either /iB or /ix is large, with deviations at smaller /i that arise because time 
scales associated with hopping rates in the same neighbourhood are not well-separated. 
We show some illustrative results in Fig. [6l the fit for the effective dynamics with 
A*B — Mt = 3 is strikingly good. However, reducing the value of /ix further has an 
effect on the diffusion front for the effective dynamics, while no effect is discernable for 
the real dynamics. This reduces the quality of the fit in this case. (The results shown 
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Figure 6. We show estimates of the function g(x) obtained by evaluating Gk{t) 
for times t within the scaling regime. Results are displayed for both real and 
effective dynamics, with /ib = 3. Symbols show results for real dynamics with 
/ix = (triangles), fi^^ = 0.5 (squares) and fiq^ = 3 (circles). The times are chosen 
to lie in the scaling regime and vary from 2^^ to 2^'*, according to the model: we 
note that while the functions g{x) are similar in all cases, the values and scalings 
of the mean square displacements are different, according to II55I I. The solid line 
shows the result for the effective dynamics and /^t = Mb = 3, while the heavy 
dashed line shows effective dynamics for fiB = 3, fiT = 0.5. The light dashed 
line is a simple exponential distribution. The scheme of DMF (not shown) gives 
reasonable agreement for fiq^ = /i^ = 3, but this agreement breaks down as )it 
(or ^b) is reduced. 



are for the physical effective dynamics scheme. For these disorder-averaged quantities, 
we note in passing that the differences between our 'physical' and 'formal' schemes 
are of the same order as the differences bewteen real and effective dynamics, with the 
physical scheme being slightly closer to the real dynamics than the formal one.) 

We also find that numerical implementation of the DMF procedure yields 
reasonable agreement with mean-square displacement and the diffusion front for the 
case fiB — — 3. However, this agreement breaks down as /ix is reduced: for 
/L^T < 1, the DMF scheme yields the wrong dynamical exponent, as discussed above. 
Thus, the main advantage of the scheme presented here is that it captures the crossover 
as gets small [see Eq. ([55])]: in this case, it also gives more accurate results for the 
propagators at fixed disorder (recall Fig. [5Jc,d,e)). 

Finally, we note that when calculating propagators, we expect the various schemes 
(DMF, Mon03, that of [10], and the one presented here) to be equivalent in the limit 
of large /i, at least at the level of the diffusion front. However, the approach to that 
limit is non-trivial and involves effects that are non-perturbative in /i: the scheme 
presented here captures some of these effects, which results in improved fits to the 
diffusion front. In the next section, we illustrate this in the case of the pure trap 
model. Of course, because of duality, an essentially identical discussion can be given 
for the pure barrier case. 

5.2. Comparison of RG schemes in the pure trap models 

We recall that the pure trap model is the case fiB 0; we then write fiT = fJ.. 
Applying the DMF method directly to the pure trap model results in a dynamical 
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Figure 7. (a) Scaling form for the disorder-averaged diffusion front for the 
pure trap model with /^t = 3, /xb = 0. We show the results for real dynamics 
and various approximations to this function, obtained by simulations of effective 
dynamics schemes. All results are taken from the scaling regime, with in particular 
t = 2^"^ for the real dynamics. We show the Mon03 scheme of Ref. 1171 . which 
yields H58| l: the scheme of [10] (labelled JSl); and the scheme discussed in this 
article (labelled JS2). (b) Distribution of barrier widths P{1) in the effective 
dynamics, scaled by its mean ^. We show results for the effective dynamics scheme 
of MonOS in which P{1) is a simple exponential, and two distributions obtained 
with the scheme JSl. 



exponent z ~ ii. As noted in [17j . this result is incorrect: the route taken by Monthus 
was simply to introduce a factor of the root mean square displacement when converting 
the rate F to a time t. With this change, the dynamical scaling given by the Mon03 
scheme [17] is correct. 

In Fig. [71 we compare the effective dynamics scheme set out in this article with 
that of [T0| and with MonOB [17]. The differences are quite striking, as we now discuss. 
In the schemes of 10^ and [17], all effective trap regions are single sites, and the RG 
scheme relates the propagator to the distribution P{1) of the widths I of effective 
barrier regions. [We normalise P{1) to dl IP{1) — 1, such that the probability for a 
randomly chosen transition state to be in a barrier region of width / is IP [I). Results 
for P{1) from our current RG scheme are not given here because even for a pure trap 
model the RG flow eventually leads to a mixture of effective traps and effective barriers 
which cannot be characterised by a single distribution P{1)] If a given barrier region 
is delimited by sites h and b + l, then the propagator is Gn,m = SnM ''+'~™ + Sn,b+i^^^^^ 
for b < m < b + l. Averaging over all initial sites m at fixed fc = n — m on a long chain 
with the relevant distribution of barrier widths, we arrive at the disorder-averaged 
diffusion front, 

G-,^ Hdl^-^Pil) (57) 
J\k\ ' 

where we have assumed that t is large, so that P(l) is smooth, and we may convert sums 
over I to integrals. Within the MonOS scheme, the distances I are all independently 
and identically distributed with an exponential form P(/) — ^~^e~'/^. This leads to 
an estimate for the diffusion front in the limit of large fi: 

« r'e-'=/« r dy—^e-y (58) 

Jo y + k/C 

so that the diffusion front is a scaling function oi x — k/a ior large times, with cr = ^. 
The assumption of Ref. [T7] is that, while working at large finite /i does affect P{1), 
these changes lead to perturbative corrections to the diffusion front. 
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However, Fig.[7]shows that the tail of the diffusion front is rather different from the 
prediction ([55]) of Ref. [T7j, at least for ^ = 3. The effective dynamics scheme discussed 
in this paper gives more accurate predictions for this tail. In fact, the convergence 
of the tail of the diffusion front to its large-/i prediction is quite slow. Instead of 
plotting the diffusion front data directly, we show barrier width distributions P{1) 
for the effective dynamics schemes of Refs. |17j and [lOj. (The distribution P{1) is 
obtained directly from the effective dynamics: in fact this was the route by which Gk 
was evaluated in Figs. [6] and [7j a).) If we define the mean barrier width to be ^, then we 
find that P{1) converges to the simple exponential distribution only if we take /i — > oo 
at a fixed value of the scaling variable Z/^. For smaller I there are corrections to this 
distribution that cannot be accounted for by treating \/ jjL perturbatively. Indeed, our 
numerical results are most consistent with the tail of the diffusion front scaling as 

logG^^ -|fc|i+" (59) 

where a > is a power that vanishes as /i — > cx). In this case it is clear that the 
limits of large ji and large k do not commute and hence that working perturbatively 
in is likely to fail when considering the large- fc limit of Gk- We have not found an 
analytical treatment which can determine the resulting exponent a, neither exactly 
nor for our effective dynamics. However, the numerical evidence of Fig. [7] is that our 
scheme does capture the non-perturbative effects which lead to slow convergence of 
the diffusion front to the large /i limit. 

6. Outlook 

In this article, we have derived duality relations that connect pairs of hopping models 
linked by an inversion of their energy landscape. The simplest case is that of models 
with equilibrium steady states and periodic boundary conditions, but we were able 
also to link models with absorbing and reflecting boundary conditions. Somewhat 
surprisingly, certain periodic systems (pure trap and barrier models with a bias) in 
which the steady state has a finite current can be analysed similarly [21]. All duality 
relations are initially expressed in terms of the relevant master operators, but we 
showed that one can then also construct the propagator of each model from its dual. 
It follows further that the disorder-averaged propagators in each pair of models are 
equal on all time and length scales. We discussed an alternative duality relation 
giving the same results, which is independent of the disorder and related to one used 
by Schiitz and Mussawisade [11] for a reaction-diffusion model. 

We have also introduced an effective dynamics scheme for these hopping models. 
It incorporates both the scheme of Ref. [T(J and that of DMF, allowing a broad class of 
models to be treated in a unified fashion. For a range of "mixed trap-barrier models" , 
including the pure barrier and trap cases, we have also shown that our scheme captures 
non-perturbative corrections to the schemes of DMF and Mon03. 

Our results also identify a few questions: can explicit expressions for the disorder- 
averaged diffusion fronts be derived for the mixed model or for the pure trap/barrier 
cases, either exactly or at least within the effective dynamics? It appears that the 
diffusion front in the mixed model depends only on the larger of /iB and ^t: can this 
be established? More speculatively, one might ask if the methods used here can be 
generalised in order to identify effective trap and barrier regions for higher-dimensional 
systems. We leave these issues for future work. 
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Appendix A. Effective dynamics 

In this appendix we give some details of the renormalisation scheme that underlies 
our effective dynamics. The scheme gives a good description of the dynamics of the 
model in the limit in which rates in the same neighbourhood are sufficiently well- 
separated. We first derive the version of the RG procedure that is most natural from 
a formal point of view, but then argue in favour of the more physically-motivated 
scheme described in the main text. Both schemes agree in the relevant limit where 
time scales are locally well-separated, and we argue that the scheme of the main text 
captures the subleading corrections to this limit more effectively. 

Appendix A.l. Formal scheme 

Suppose that we have a master operator T4^R(r) of the form given in ([50]) . and that 
this operator gives an accurate description of motion on time scales longer than 
r~^. We wish to construct a projection operator P{V — 5V) which represents a 
good approximation to the operator Pc^{T — 5T) of (|T7)) . so that M^R,(r — 5T) = 
P(r — i5r)M^R(r)P(r — SV) gives an accurate description of motion on time scales 
longer than (F — (JF)^^. In addition, for the scheme to represent a renormalisation 
group flow, we require that Wr(F — ST) is also of the form given in ([50]) . 

As discussed in the text, we begin by estimating an eigenmodc of lyR(F) that is 
concerned with fast motion. To achieve this, we imagine that all the rates r^ and 
are associated with very slow motion, except the triplet {r^, In that case 
we can write the master operator as 

/ ^„ \ / \ 

WR(r) « Wo = (|P„-i) -(C+r„) (Q„| (A.l) 

Diagonalising yields three eigenvalues. Since we have assumed that transitions 
out of site a — 1 are very slow the 'steady state', i.e. the right eigenvector with 
eigenvalue zero, of this reduced system is simply localised on that site. Then 
there are two negative eigenvalues whose moduli are p± = \{ra + + ^a+i) ± 
^y/fa + 2ra(£a + ^a+i) + {(-a — ^a+iY ■ These eigenvalues are associated with fast 
(-f) and slow (— ) motion. We identify p+ as a rate for fast motion to the right 
from effective trap a. The effective dynamics proceeds by successive removal of the 
fastest such modes: recall where p+ is written as pa- As discussed in the main 
text in addition to triplets of rates (r^, ^a, £q,+i), we also consider triplets such as 
(£q, Ta, ra_i), for which the same treatment applies, with the rate Aq of the fastest 
mode given in ([39| . In the discussion below we assume for concreteness that the 
largest approximate eigenvalue among the {/Cq, Aq,} is Pa = p+- 

To accomplish the removal of the fastest mode, we will project the original master 
operator Wr(F) onto the basis spanned by the eigenvectors associated with the slow 
motion. More precisely, the zero eigenvectors of the approximate master operator Wq 
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are \Pa-i) to the right and (esl — (Qq-iI + {Qa\ + {Qa+i\ to the left. The right and 
left slow eigenvectors, corresponding to eigenvalue /9_, we write as \p-) and (/0_|. We 
therefore define a projection operator whose matrix elements will give the propagator 
on time scales longer than 

a — 2 oo 

V+= \Pc.'){Qc.'\ + \Po.-i){e;\ + \p-){p-\+ E l^"')(Qa'| (A.2) 

a' — — oo a' — a+2 

In general, the operator V+Wyi{V)V+ is not of the same form as WR(r): it 
contains hopping processes between next-nearest neighbours for the effective traps, 
and does not represent a suitable approximation to WR(r — (5r). This effect is familiar 
in renormalisation schemes, and requires irrelevant terms in the master equation to 
be discarded. In our situation, next-nearest neighbour hopping becomes irrelevant 
as time scales become well-separated. For example, in the barrier-like case where 
fa ~ £a+i but ^ £a, we have 

a-2 

V+^V^= \Po.'){Qo.'\ + \Pc.-i){Qc.-i\ + \P'){{Q^\ + {Qc.+i\) 

a' — — oo 
oo 

+ (A-3) 

q'=q+2 

with 

We identify \P') as the right eigenvector associated with a new effective trap that 
combines traps a and a -\- 1. The resulting master operator VBW^iT)!^^ is now of 
the same form as WR(r), with the free energy of the new effective trap given by 
e^^.a+i = e-*^" -I- e^"*^ . Thus, we may perform an RG step by taking V{V — 5T) = Pb 
as our approximation to "Poxlr — 5T). Such an RG step is illustrated by the transition 
between Fig.[5}3 and Fig. [31;: consistent with that figure, the transition states on either 
side of the new trap are unchanged during this procedure. 

On the other hand, in the trap- like case where (.^ ~ r^. but ^ ^a+l^ we have 

P+«Pt= Y l^"')(Q"'l + l^o-i)((e3|-(Q'|) + |Pa+i)(Q'| 



E 



|Pa')(Qa'l (A.5) 



with 



= — ^(Qa| + (ga + l|. (A.6) 

As for the previous case, 7^TWR(r)'PT is of the same form as WR(r), so the choice 
■p(r — 5T) = Vt corresponds to a valid RG step. In this step, trap a has been 
incorporated into a new effective barrier region that merges the old barrier regions 
a — i and a + ^ • Its escape properties to the remaining effective traps a+1 and a — 1 
are described by the eigenvectors {Q'\ and (63] — {Q'\ — {Qa-i\ +ia/{'''a + ^a){Qa\: 
respectively. An example of such an RG step takes place between Figs. [3^ and[3)D. 

It should be noted that if Tq. is much greater than both ia and £a+i then 
both (|A.3P and (jA.Sp reduce to the case of DMF, which is \P') = \Pa+i) and 
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{Q'\ = {Qa\ + {Qa+i\- In practice, for any renormalisation step, we choose either 
to combine traps by taking P(r — 5T) — or otherwise to combine barriers by 
taking P{T — ST) = Px- In either case the resuhing VFr (F — ST) is indeed of the same 
form as VFFi(r): this ensures that our procedure is a vahd renormahsation flow in the 
space of hopping models. 

As discussed in Sec. 14. 4| the vahdity of the renormahsation scheme requires that 
the projection operator evolves with F such that V{T — ST) ss Pe^{T — ST). Assuming 
that we combine traps a and a + 1, we should have Pb ~ P+ ~ Pcx- Applying 
perturbation theory to the fast eigenvectors, we find that corrections are small if 
fa + ia+i ^ ra-ijlaTfa+i, £a+2- Similarly, if we combine barriers a ± |, we require 
^a + fa ^ ia-i^fa-i, ia+i,ra+i- Essentially, if the rate Tq, is large at a given stage of 
the RG (in the sense that pa > T — ST), then must be larger than all rates in the 
neighbourhood, except for either ia or ia+i, one of which may be comparable to ra- 

Appendix A. 2. Illustrative example, and comparison with physical scheme of the 
main text 

The formal scheme derived above differs from the one set out in the main text in the 
way effective traps and barriers are combined. We illustrate this with the system of 
Fig. [H The bare master operator is 

W = --- + m~ - ro(0|) + (|1) - |2))(^2(2| - ri(l|) + • • • (A.7) 

where the • ■ • indicate the remaining terms in the master operator, including those for 
hopping into and out of this segment of the chain. 

Applying the first step of the formal RG scheme, we construct the projection 
operator 

P = --- + |0)(Qo| + |2)(g2| + --- (A.8) 

E 1 ~F E 3 —F 

where (Qol = (0| +e * (1| and (Q2I = (2| + e « (1| are constructed according to 

(|A.6p , with e-*^ = e i + e i . This leads to the same result as the physical scheme 
of the text, since (IA.8|) is also consistent with ([53]) above. The renormalised master 
operator is 

= PWP = . . . + (|0) - |2))e-^(e-^^(g2| - e-^°(Qo|) + • • • (A.9) 
Applying the formal scheme again to this operator, the projection operator is 

v'= •■• + |Po2)((goi + (Q2i) + --- 

= ••• + I^02>((0| + (1| + (2|) + --- (A.IO) 

where IP02) = e^^^^ (e^" |0) + e^'^\2)) with e-'^"^ = e-^" + e-^^ According to the scheme 
of the main text, we would have obtained a similar result, but with the replacement 
\Po2) \P012) =e-^°i2(e^''|0) + e^i|l) + e^=|2)), where e^°i= = e^» + e^^ + e^^ 

In the limit where rates are well-separated and the formal scheme is is exact, it 
follows from that Ei <C -Eq, -E'2i and in this case, the physical and formal schemes 
coincide. Indeed, it may be shown that the errors associated with the physical and 
formal schemes are of the same order. We conclude that the physical scheme is at 
least as appropriate as the formal one and indeed it can be verified that the errors are 
of the same order. 

The key point is that the physical scheme of the text is based on the assumption 
of equilibration within effective traps. Under that assumption, it is clear that the 



Duality symmetries and effective dynamics in disordered hopping models 28 

states \Pa) should coincide with Boltzmann distributions over the sites within the 
trap. In general, the formal scheme gives a state \Pa) that is finite only on a subset 
of the sites within the trap. This is a necessary feature of the formal scheme, because 
the \Pa) at each stage are constructed from the \Pa) of the previous stage. In the 
language of section |4?T1 the slow co-ordinates (i) must be linear combinations of the 
Pa of the previous stage. Thus, since site 1 in the example of Fig. [3] is not contained 
in any of the pa after the first step of the RG, it can never be part of an effective 
trap at any future stage. Equivalently, the \Pa) can have no contribution from site 
1 in future stages, and may not correspond to Boltzmann distributions of the whole 
trap. This in turn means that the approximate propagator (n|7^(r)|m) is equal to 
zero for some final sites n within effective traps, underestimating the true value of 
G„m(r~^) there. The differences between the formal and physical schemes are shown 
numerically in Fig. [5l which illustrates how the formal scheme underestimates the 
probability of propagation onto certain sites within the effective traps, for specific 
disorder realisations. As discussed in Sec. [51 the physical and formal schemes show 
differences in the tails of the disorder-averaged diffusion front that are of the same 
order as the deviations between real and effective dynamics. 

Finally, we note that the physical scheme is not strictly a renormalisation flow in 
that the rates and free energies at a given stage, depend not just on the renormalised 
operator Wr at the previous stage, but on all the bare energies Ei and E^_^i. If 
this feature is considered undesirable, the formal scheme may be used. However, the 
interpretation of the physical scheme as a partition into contiguous trap and barrier 
regions and the intuitive idea of equilibration within effective trap regions means that 
we prefer that route. 
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